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The extension of thermodynamics into the quantum regime has received much attention in recent 
years. A primary objective of current research is to find thermodynamic tasks which can be enhanced 
by quantum mechanical effects. With this goal in mind, we explore the finite-time dynamics of 
absorption refrigerators composed of three quantum bits (qubits). The aim of this finite-time cooling 
is to reach low temperatures as fast as possible and subsequently extract the cold particle to exploit 
it for information processing purposes. We show that the coherent oscillations inherent to quantum 
dynamics can be harnessed to reach temperatures that are colder than the steady state in orders of 
magnitude less time, thereby providing a fast source of low-entropy qubits. This effect demonstrates 
that quantum thermal machines can surpass classical ones, reminiscent of quantum advantages in 
other fields, and is applicable to a broad range of technologically important scenarios. 


I. INTRODUCTION 

The development of classical thermodynamics in the 
19th century underpinned the Industrial Revolution, and 
the enormous economic growth and social changes that 
followed. Now, in the 21st century, the burgeoning quan¬ 
tum technological revolution promises unprecedented ad¬ 
vances in our computation and communication capabili¬ 
ties, enabled by harnessing quantum coherence. As our 
machines are scaled down into the quantum regime, it 
is of prime importance to understand how quantum me¬ 
chanics affects the operation of these devices. This prob¬ 
lem has attracted great interest to the field of quantum 
thermodynamics over the last few years. 

One useful approach in this regard is to explore simple 
physical models which highlight novel aspects of quantum 
thermal machines. The quantum absorption refrigerator 
is the quantum extension of a classical machine devised in 
the 19th century (see, for example. Ref. [1] and references 
therein). The smallest possible model with couplings be¬ 
tween physical particles and thermal reservoirs was first 
studied by Linden et al. [2], who considered a three- 
qubit refrigerator. While the history of studying these 
machines dates back a long time even in the quantum 
regime, this three-qubit model was the first where the 
role of quantum information resources was studied [3], 
revealing that entanglement in the steady state prohibits 
achieving perfect Carnot efficiency, but potentially in¬ 
creases cooling efficiency. A wide variety of other quan¬ 
tum absorption refrigerator models have also been pro¬ 
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posed in the recent literature [4-8]. 

Designing thermodynamic processes that can be en¬ 
hanced by quantum dynamics is a pivotal challenge in the 
field of quantum thermodynamics. One of the paradig¬ 
matic thermodynamic tasks concerns work efficiency at 
the quantum scale. Here already the very definition of 
quantum mechanical work is debated [9-11], yet in differ¬ 
ent scenarios quantum mechanical advantages seem pos¬ 
sible [12-15]. Another avenue in this endeavour is the 
engineering of the environment itself to enhance quan¬ 
tum processes [16-19]. We circumvent the potential con¬ 
troversies regarding the practical value of work and ef¬ 
ficiency by concentrating on a different figure of merit, 
namely the achievable temperature, and by using ther¬ 
mal baths to drive the refrigerator, therefore needing no 
notion of work. 

Previous research on the quantum aspects of heat en¬ 
gines and refrigerators has focused almost exclusively 
on their operation in the steady state. However, in 
many applications one wishes to reach low temperatures 
as rapidly as possible, in which case understanding the 
short-time behaviour is essential. In particular, the cool¬ 
ing may be applied only transiently, after which the cold 
object is extracted for use. As a somewhat frivolous yet 
illustrative example, consider the problem of refrigerat¬ 
ing a beverage on a hot day. Maximum enjoyment is ob¬ 
tained if the beverage can be cooled and then consumed 
quickly, at a temperature significantly lower than that of 
the environment. A more serious example could be the 
initialisation of a register of qubits for quantum informa¬ 
tion processing, where the aim is to produce states with 
high purity (low entropy). Fast cooling is advantageous 
here since it may reduce the overall time taken to com¬ 
plete the quantum information protocol. Both of these 
situations exemplify what we call single-shot cooling: the 
one-time application of a refrigeration device in order to 
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considerably and rapidly reduce the temperature of the 
object in question. 

In the present work, we study the application of three- 
qubit absorption refrigerators to single-shot cooling. By 
considering variations of the basic processes underlying 
energy dissipation and transport, we elucidate the role 
of coherence in the operation of such a device. Further¬ 
more, we demonstrate that dramatic improvements can 
be obtained in both the cooling time and the achievable 
temperatures by taking advantage of coherent oscillations 
that appear in the transient dynamics of the refrigerator. 
Most importantly, the presence of coherence in the device 
allows one to reach lower temperatures than the steady 
state. This use of quantum coherence in order to cool be¬ 
low the steady-state temperature we refer to as quantum 
single-shot cooling. 

Note that several physical implementations of quan¬ 
tum thermal machines have recently been proposed [20- 
25 ], and experimental efforts to construct such devices 
in the laboratory are currently under way. Our find¬ 
ings could be relevant for future experiments, which may 
be practically limited by finite coherence times. From 
a more fundamental perspective, our scheme provides 
one of the first examples in which quantum coherence 
plays an active and necessary role in improving the per¬ 
formance of machines driven only by thermal noise. 



FIG. 1 . Cartoon illustrating our theoretical set-up. (a) A 
stream of thermal qubits passes sequentially through a quan¬ 
tum absorption refrigerator. Each thermal qubit interacts 
with the refrigerator for a finite time, and exits the device 
at a lower temperature, (b) The canonical model of ther¬ 
mal reservoirs (Model I) consists of an infinite collection of 
harmonic oscillators, characterised by a bandwidth Q and a 
dimensionless coupling strength a. (c) An alternative bath 
model (Model II) consists of a single fictitious qubit undergo¬ 
ing Markovian gain and decay processes, characterised by a 
linewidth 7 and a dimensionless coupling strength 77. 


II. THEORETICAL FRAMEWORK 
A. Description of the refrigerator 

We consider a quantum absorption refrigerator com¬ 
prising three qubits described by standard Pauli opera¬ 
tors with the local Hamiltonian 

1 ^ 

Uoc = -Y.Eial. ( 1 ) 

i=l 

Throughout this article, we employ units of energy, time 
and temperature such that Ei = h = 1 and /cb = 1- 
The qubits are coupled together according to a three- 
body interaction 

H = ^/|010)(101|+h.c., (2) 

where the computational basis states |0), |1) denote the 
eigenvectors of In order for this interaction to be 
energy-conserving, in the sense that [i^ioc, H] = 0, we de¬ 
mand that the qubit energies satisfy E2 = Ei E^. The 
interaction (2) then drives resonant transitions within the 
transport subspace spanned by the states |010) and |101). 

Each qubit i is in contact with an independent heat 
bath which drives it towards a thermal equilibrium 
state at inverse temperature jSi = 1 /T^, where Ti < T2 < 
T3. When the system parameters are judiciously chosen, 
the dynamics propels qubit 1 (the “cold qubit”) towards 
a new temperature Ti < Ti. It is in this sense that the 


system behaves as a refrigerator, whose basic operating 
principle may be understood as follows. The interac¬ 
tion Hamiltonian ( 2 ) couples the cold qubit to a “virtual 
qubit” [ 26 ] comprising the {|01),|10)} subspace of the 
Hilbert space of qubits 2 and 3 . When these qubits are 
in equilibrium with their respective baths, the popula¬ 
tions of the virtual qubit states are thermally distributed 
at an effective temperature given by 


E2 — Es 

(32E2 — l^sEs 


( 3 ) 


If the parameters of the refrigerator are arranged so that 
Ty < Ti, then the cold qubit equilibrates to a lower tem¬ 
perature as it exchanges energy with the virtual qubit on 
the approach to the steady state. 

In the following, we specialise to the case where Ti = 
T2 =: and we define := T3. This is in many ways 

the most natural scenario, where we have a single source 
of free energy — namely the hot bath at temperature 
Th — enabling us to cool below the ambient “room tem¬ 
perature” Ty. So long as qubits 1 and 2 are sufficiently 
spatially separated [ 27 ], we can treat the effect of their 
common environment as arising from two independent 
baths [ 28 ]. 

The general set-up that we have in mind is illustrated 
in Fig. 1 (a). We assume that all three qubits are initially 
in thermal equilibrium with their respective reservoirs. 
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The initial state is therefore 

Pa{0) = -JY ’ 


(4) 


where pA {t) denotes the reduced state of the refrigerator 
qubits at time t. At t = 0 the interaction ( 2 ) between 
the qubits is switched on. The aim is then to prepare the 
cold qubit in a low-entropy state. (We emphasise that 
this is quite different from the usual objective consid¬ 
ered by previous authors, namely to maintain the cold 
qubit at a low temperature.) Once the lowest achievable 
temperature has been reached at time to, our cooling pro¬ 
tocol is complete, and the cold qubit may be extracted 
for use, for example, in quantum information processing. 
Classically, we usually expect that the lowest achievable 
temperature occurs only as t ^ oo. The existence of 
a finite optimal extraction time to is the essence of the 
quantum single-shot cooling protocol. If we suppose now 
that there exists a large supply of qubits thermalised at 
the ambient temperature T^, we can perform the same 
procedure repeatedly to produce a steady stream of low- 
entropy qubits. 


B. Models of thermalisation 

In order to quantitatively analyse the refrigerator we 
must specify a thermalisation model. On the other hand, 
one would like to obtain general results that are indepen¬ 
dent of any particular model. In order to avoid being too 
restricted by our assumptions, we employ two different 
approaches to modelling the thermal baths. In each case, 
we assume that the baths are Markovian (memory-less) 
and weakly coupled to the refrigerator, and we describe 
the dynamics by a master equation in Lindblad form. In 
the following paragraphs, we describe these models in a 
qualitative way, deferring the full details to Appendix A. 


tained when Q is much larger than all other frequency 
scales and a ^ 1 , so that the dissipation rates are much 
smaller than the natural frequencies {Ei} of the qubits. 

Some care must be taken when treating the effect of 
the inter-qubit coupling on the thermalisation dynamics. 
We are able to derive two different master equations de¬ 
pending on the magnitude of g. The first equation is valid 
in the strong-coupling limit, where g is much larger than 
the dissipation rates [29]. The second master equation 
holds in the weak-coupling limit, when g is comparable 
to or smaller than the dissipation rates [30, 31]. 


2. Model II 

An alternative bath model consists of assigning to each 
qubit of the refrigerator an additional fictitious qubit, 
which is damped by a perfectly Markovian thermal bath 
at temperature T^, such that the spontaneous emission 
rate is 7 (Fig. 1(c)). Excitations are then allowed to hop 
between each fictitious qubit and its associated physical 
qubit at a rate 777 . This simulates an effective thermal 
reservoir with a memory time of order 7 “^, coupled with 
strength 7 to the physical qubit. Markovian dynamics of 
the physical qubits are therefore obtained when 7^1 
and 7 is larger than all other frequency scales. The en¬ 
ergy splitting of each fictitious qubit is chosen to be res¬ 
onant with the frequency Ei of the physical qubit, which 
ensures that each physical qubit in isolation is driven to¬ 
wards a local thermal state at temperature T^. 

This model enables us to explore the strong-coupling 
limit of large g without invoking the rotating-wave ap¬ 
proximation, which involves a time-averaging over the 
autonomous dynamics of the refrigerator. This time¬ 
averaging has been shown to lead to unphysical predic¬ 
tions in certain scenarios where the open system is cou¬ 
pled to multiple heat baths at different temperatures [30] . 
Thus Model II provides us with an independent check on 
the validity of our results. 


1. Model I 


III. RESULTS 


The canonical procedure to model Markovian reser¬ 
voirs is to suppose that each qubit is coupled to an infi¬ 
nite collection of harmonic oscillators spanning a broad 
range of frequencies (Fig. 1(b)). In our case, the baths 
are described by identical Ohmic spectral functions of the 
form 


J(ijj) = auje (5) 

This function quantifies the strength of the coupling be¬ 
tween each qubit and the oscillators near frequency cj, 
weighted by the density of states of the reservoir (see 
Appendix A). The effect of the baths is therefore char¬ 
acterised by two parameters: a dimensionless coupling 
strength a, and a frequency cut-off U leading to a bath 
memory time of order Markovian dynamics is ob- 


A. Short-time dynamics of the refrigerator 

We now present our quantitative results, obtained by 
numerical solution of the equations of motion. Our first 
observation is that sufficiently strong coherent coupling 
between the qubits drives damped Rabi oscillations of 
the local qubit populations. However, unlike Rabi oscil¬ 
lations due to local driving fields, the three-body inter¬ 
action ( 2 ) does not induce any local coherences between 
the qubit populations. The reduced state of each qubit 
is diagonal at all times and may therefore be assigned an 
effective temperature 


fS) 


-Ei 

2tanh“^((Tf (t)) 


(6) 
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To illustrate this feature of the short-time dynamics, 
we plot several examples of the evolution of the cold 
qubit temperature in Fig. 2, which demonstrate that the 
Rabi oscillations allow the cold qubit to reach lower tem¬ 
peratures than the steady state. In the strong-coupling 
regime we find damped temperature oscillations with the 
approximate period tt/^. Therefore, the optimal quan¬ 
tum single-shot cooling procedure consists of extracting 
the cold qubit after a time to ~ 7r/(2^). So long as the 
coupling g is larger than the relaxation rate, the same 
qualitative behaviour is found in both Model I (solid line 
in Fig. 2(a)) and Model II (solid line in Fig. 2(c)). In 
contrast, when g is much smaller than the relaxation 
rate, these oscillations are over-damped so that no tem¬ 
perature minimum occurs in a finite time (solid line in 
Fig. 2(b)). Nevertheless, we have found that quantum 
single-shot cooling is possible over a very broad range 
of parameters, so long as the coupling g is significantly 
larger than the relaxation rate. 

The physical origin of the temperature oscillations is 
the exchange of energy between the refrigerator qubits 
due to the interaction (2). Our second important obser¬ 
vation is that this energy transport is driven by coherence 
in the transport subspace. This can be seen straightfor¬ 
wardly by examining the Heisenberg equations of motion 
for the local energy expectation values hi = Ei{a^)l2. 
The resulting expression for the cold qubit is of the form 

^ = Qiit)-2gE,lmC{t). (7) 

This equation represents an energy balance between the 
rate of heat absorbed from the bath Qi (t) and the co¬ 
herent flow of energy into the other qubits, which is pro¬ 
portional to the imaginary part of the coherence in the 
transport subspace: 

C(t)=lV[p(i)|010)(101|], (8) 

where p{t) denotes the quantum state. Eqs. (7) & (8) give 
a direct link between the presence of coherence and the 
flow of energy across the refrigerator. We expect an anal¬ 
ogous relationship to hold for any quantum refrigerator 
in which energy transport between distinct sub-systems 
occurs via coherent Hamiltonian evolution. 

Note that the same mechanism, whereby coherence en¬ 
hances the flow of energy, has recently been shown to lead 
to an operational advantage for heat engines and power 
refrigerators working in the steady state [15]. We also 
mention the similarity of our protocol with algorithmic 
cooling [32, 33]. Indeed, the three-qubit unitary swap op¬ 
eration we employ is formally identical to that proposed 
for specific algorithmic cooling protocols [34] , which may 
also be used in a single-shot operation similar to our pro¬ 
posal. The novel feature of our set-up is that no work 
is required to perform cooling. Rather, the free energy 
source is provided by the temperature difference between 
the hot thermal bath at temperature and the ambient 
environment at temperature T^. This temperature dif- 




t t 
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FIG. 2 . Effective temperature dynamics of the cold qubit, 
with El = 1, E 2 = 2, Tr = Ti = T 2 = 50 and = T 3 = 100. 
(a-c) Blue solid lines indicate the evolution of an initial ther¬ 
mal product state, red dotted lines depict the evolution with 
1 % of the maximum coherence added to the initial state, 
(a) Model I in the strong-coupling regime, (b) Model I in 
the weak-coupling regime, plotted with the corresponding 
stochastic model (black dot-dashed line), (c) Model II in the 
strong-coupling regime, (d) Short-time effective temperature 
dynamics of an ensemble of cold qubits, where coherence with 
some phase noise is added to the initial state. The magnitude 
of the initial coherence is 5% of the maximum and the phase 
uncertainties are (50 = 0 (black solid line), (50 = 7 r /4 (blue 
dotted line) and (50 = tt (red dot-dashed line). The tempera¬ 
ture minimum in the absence of phase noise is shown by the 
green vertical dashed line. The data shown in (d) are calcu¬ 
lated using Model I with the same parameters as (a). 


ference establishes a population bias between the states 
1101) and |010), which enables the unitary generated by 
V to redistribute excitations, thereby cooling the target 
qubit. 


B. Effect of initial coherence 

To further elucidate the fundamental role of coherence 
in energy transport, we consider a situation where some 
coherence in the transport subspace is added to the initial 
state, without modifying the thermally distributed popu¬ 
lations. In order to ensure the positivity of the quantum 
state, the magnitude of this coherence is upper bounded 
by |C(0)| < Cmax, where 

Cmax = n ^sech • (9) 

Eor each example in Eigs. 2 (a-c) we have also plotted 
the dynamics with a very small amount of initial coher¬ 
ence C(0) = iCmax/100- We find that the amplitude of 
the temperature oscillations is noticeably enhanced in 
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all three cases. In the weak-coupling example plotted 
(red dotted line in Fig. 2(b)), the initial coherence gives 
the most dramatic advantage, since now the tempera¬ 
ture minimum occurs in finite time. Initial coherence is 
actually a necessary ingredient for quantum single-shot 
cooling in this case and thus should be considered as an 
additional non-classical resource to the free energy of the 
hot bath. 

This result reinforces the notion that coherence is a 
useful resource in quantum thermodynamics, which has 
lately been highlighted by a number of other authors. 
Indeed, it has recently been shown that coherence can 
improve the performance of thermal machines in the al¬ 
ternative context of a time-dependent Hamiltonian [35]. 
In the abstract resource-theoretic approach to quantum 
thermodynamics, coherence is a resource which is non¬ 
increasing under any allowed operation [36, 37], and from 
which work can be extracted [13, 38]. Coherence may 
also play the role of a catalyst, but only for infinite¬ 
dimensional quantum systems [39]. Furthermore, it has 
been demonstrated explicitly that work can be extracted 
from coherence by quantum measurements [40]. 

Note that, according to Eq. (7), the phase of C(0) de¬ 
termines the direction of the initial flow of energy into the 
qubit. Therefore, it is even possible, for example, to heat 
the qubit by adding initial coherence with the opposite 
phase, even when the system behaves as a refrigerator 
in the steady state. Likewise, one can transiently cool 
even when the steady-state behaviour is that of a heat 
pump. Similar phase effects have recently been described 
in Ref. [41], where the authors show that bath fluctua¬ 
tions can revert the detailed balance condition in certain 
open quantum systems, creating a net flux of energy from 
the environment into the system. 

Due to the aforementioned sensitivity of the dynamics 
to the phase of the coherence, it is important to con¬ 
sider the impact of any phase uncertainty resulting from 
imperfect initial state preparation. Specifically, we sup¬ 
pose that the initial coherence is given by C(0) = ire^"^, 
where 0 is a zero-mean random variable. Assuming that 
the cooling protocol is performed repeatedly on multiple 
qubits, the relevant quantity to consider is the tempera¬ 
ture of the ensemble of cold qubits. For the rest of this 
section it should be understood that the term tempera¬ 
ture refers to the property of the ensemble. (A charac¬ 
terisation of the temperature fluctuations of individual 
qubits is beyond the scope of this article.) 

The introduction of phase noise leads to two effects: an 
increase of the minimum temperature, and a shift in the 
time at which this minimum occurs. In Appendix B we 
provide an approximate analytical quantification of the 
impact of phase noise, assuming that the phase fluctu¬ 
ations are small and approximately Gaussian. We have 
also numerically investigated the case of non-Gaussian 
phase noise, in particular the case where (j) is uniformly 
distributed in the range (j) G [—(50/2,-h(^0/2]. Some ex¬ 
amples of the resulting temperature dynamics of the cold 
qubit ensemble are plotted in Fig. 2(d). 


We find that a narrow phase noise distribution leads 
to a small shift in the minimum temperature (e.g. blue 
dotted line in Fig. 2(d)). The most deleterious impact is 
obtained when the phase noise distribution is very broad, 
so that the temperature minimum is shifted to a much 
later time. Extracting the cold qubit at the expected 
minimum may actually lead to heating of the cold qubit 
(red dot-dashed line in Eig. 2(d)). On the other hand, 
it may be possible to predict the dynamics of the en¬ 
semble temperature in advance, e.g. because the phase 
noise distribution is known. In this case it is possible 
to extract the cold qubit at the appropriate time to ob¬ 
tain some temperature advantage, unless this time is so 
late that thermal dissipation damps away the coherent 
oscillations. 


C. Stochastic absorption refrigerator 

In the previous two sections we presented strong evi¬ 
dence that coherence is a useful resource that may be har¬ 
nessed to produce a significant advantage to single-shot 
cooling. On the other hand, we can show that coherence 
is not a necessary ingredient for quantum absorption re¬ 
frigerators operating in the steady state. This is because 
energy transport can also be described by a stochastic 
process in which excitations are transferred incoherently 
between the qubits. Such a classical description is ap¬ 
propriate in the presence of strong dephasing [42] , which 
in our case comes directly from the thermal baths (see 
Ref. [15] for an alternative dephasing model in the con¬ 
text of power refrigerators). 

As a concrete example, in Appendix G we derive an 
effective master equation describing the asymptotic re¬ 
laxation of the populations in the computational basis, 
valid when the coupling g is much smaller than the dissi¬ 
pation rates. We show that energy transport in this limit 
can be effectively modelled by stochastic transitions be¬ 
tween the states |101) and |010) at a rate 2^^/r, where 
F is a characteristic dephasing rate due to the action of 
the thermal baths (see Eq. (G14)). 

The effective master equation accurately describes the 
long-time temperature dynamics of the cold qubit on its 
approach to the steady state (black dot-dashed line in 
Eig. 2(b)). This result demonstrates that three-qubit re¬ 
frigerators are capable of cooling even when their dynam¬ 
ics is entirely classical. We see from Eig. 2(b) that there 
is no difference between the quantum and classical model 
in terms of the temperatures which can be achieved in the 
steady state. 

In the transient regime, on the other hand, the temper¬ 
ature dynamics are qualitatively different in the quantum 
and classical cases. Glassical refrigeration models de¬ 
scribed by stochastic rate equations generically exhibit 
pure exponential relaxation, precluding the possibility of 
single-shot cooling below the steady-state temperature. 
This lies in stark contrast to the generic behaviour we 
have found in strongly coupled quantum absorption re- 
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FIG. 3. (a) Example evolution using the same parameters 
as Fig. 2(a), with (blue solid line) and without (green dot¬ 
ted line) a switch-off of the interaction at time to = 7r/(2^), 
illustrating the quantum temperature and time advantages 
AT and tg, respectively. The steady-state temperature Too 
is shown by the red dashed line, (b) Trade-off between the 
fractional temperature advantage and the time advantage (in 
units of the relaxation time yF^)? with the same parameters 
as (a) apart from the temperatures. Each line shows the re¬ 
sults for a fixed room temperature Tr and a range of hot bath 
temperatures Th G [Tr^H-l, Tr^ + 200]. The value of tq decreases 
as Th is increased in each case. 

frigerators, i.e. temperature oscillations enabling us to 
rapidly achieve temperatures lower than the steady state. 

D. Comparison with steady-state cooling 

Returning now to the fully quantum case, we would 
like to compare the performance of quantum single-shot 
cooling and steady-state refrigeration. Let us define the 
temperature advantage AT = To© — Ti{to) as the dif¬ 
ference between the steady-state temperature To© of the 
cold qubit and the quantum single-shot cooling minimum 
Ti(to). This is an important figure of merit quantifying 
the advantage gained from quantum single-shot cooling 
compared to steady-state cooling. Once the cold qubit 
is decoupled from the refrigerator at time to, its temper¬ 
ature increases as it equilibrates with the environment. 
The advantage gained from quantum single-shot cooling 
lasts only until the temperature of the cold qubit grows 
larger than the steady-state temperature, which occurs 
at a time ti. This motivates us to define the quantum 
advantage time tq = ti — to, which is another impor¬ 
tant quantity characterising the performance of the quan¬ 
tum single-shot refrigerator compared to its steady-state 
counterpart. See Fig. 3(a) for a graphical depiction of 
these quantities. 

For simplicity, in this section we work in the strong¬ 
coupling limit of Model I, where there always exists some 
quantum single-shot cooling advantage and the optimal 
switch-off time is given approximately by to ~ 7r/(2^). 
In this case, the quantum advantage time tq tends to 
decrease as either of the bath temperatures Tr and Th 
is increased. Of course, as the room temperature is var¬ 
ied, the rate of equilibration of the cold qubit with its 
environment changes, which in turn affects the value of 
tq. In order to remove this trivial dependence on T^, 
we measure the quantum advantage time in units of the 


relaxation time 7 “^, where 7 ^ = [ 71 (^ 1 ) + 7 i(—Ti )]/2 is 
the arithmetic mean of the thermal gain and decay rates 
acting on the cold qubit in isolation [43] (the explicit 
expressions for 7 i(±Ti) are given in Appendix A). 

For a given set of fridge energies {Ei^g} and a fixed 
room temperature T^, one would like to optimise the tem¬ 
perature of the free energy source Th to obtain the largest 
fractional temperature advantage AT/T^o and also the 
longest possible advantage time tq. We show in Fig. 3(b) 
that there exists a trade-off between these two quanti¬ 
ties, which cannot be simultaneously maximised. This 
trade-off is reminiscent of the competing goals of effi¬ 
ciency and power when optimising thermal machines that 
operate in the steady-state. Interestingly, there exists a 
“sweet spot” where the fractional temperature advan¬ 
tage is maximised, while the advantage time remains ad¬ 
equately large. This immediately suggests a favorable op¬ 
erating regime in which the quantum single-shot cooling 
protocol is preferable to steady-state cooling. However, 
in a practical setting the optimum refrigerator param¬ 
eters will of course depend on the specific application 
in question. Finally, it is worth noting that introduc¬ 
ing coherence in the transport subspace into the initial 
quantum state further enhances the fractional tempera¬ 
ture advantage for a given tg, thus coherence is also a 
useful resource in this context. 


IV. CONCLUSIONS 

We have studied the dynamical evolution of three- 
qubit absorption refrigerators in different regimes, using 
two different models of thermal dissipation in open quan¬ 
tum systems. In both these models we encounter oscil¬ 
lations in the temperature of the target qubit below the 
steady-state temperature, reinforcing the notion of a uni¬ 
versal and robust feature of quantum refrigerators. We 
show that with the right timing these oscillations can be 
exploited in a quantum single-shot cooling protocol to 
yield a constant stream of cold qubits. Numerous quan¬ 
tum information processing protocols require (approxi¬ 
mately) pure input states [44], potentially making quan¬ 
tum thermal machines a useful addition as pre-cursors to 
information processing protocols. 

Experimental implementation of quantum single-shot 
cooling calls for the strong-coupling limit, in which the 
rate of coherent energy transport exceeds the rate of ther¬ 
mal dissipation. This parameter regime may be difficult 
to achieve in currently proposed set-ups [21-23]. Never¬ 
theless, we hope that our results will serve as an encour¬ 
agement for further theoretical and experimental work 
towards achieving strong coupling in quantum absorp¬ 
tion refrigerators. Note that a recent work has found 
evidence that quantum single-shot cooling is even possi¬ 
ble in the weak-coupling limit [45], providing yet more 
motivation to experimentally study transient effects in 
quantum thermal machines. 

The advantage of the absorption refrigerator lies in its 
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high degree of autonomy: no external energy source is 
required to keep the machine running. Since the initial 
states we consider come for free in the context of thermo¬ 
dynamic resource theories [46], or more explicitly come 
directly from their respective baths, two thermal reser¬ 
voirs at a given temperature are sufficient to let the ma¬ 
chine run for as long as their temperatures do not change 
significantly (which for macroscopic baths in contact with 
quantum systems is sufficiently long for all practical pur¬ 
poses). 

Note that our scheme is not completely autonomous, 
because we have assumed that the qubit-qubit interaction 
can be switched on and off by external control. However, 
it is also possible to envisage a scenario where the switch¬ 
ing is performed by a quantum clock [47]. In this case, 
the global Hamiltonian would be time-independent and 
the protocol becomes fully autonomous. However, a more 
detailed study is required in order to understand how a 
quantum mechanical clock may be used to control the 
machine: this will form the topic of a future publication. 

Apart from the prospect of information processing, our 
results help to elucidate the quantum nature of thermo¬ 
dynamics. Compared to its classical counterpart, ther¬ 
modynamics at the quantum scale often results in ad¬ 
ditional constraints and limitations [10, 48] due to the 
discretised nature of the fundamental states. The cen¬ 
tral question in this context concerns the actual impact 
of coherence, entanglement and other genuine quantum 
features on the potential transformations in quantum 
thermodynamics. Indeed, researchers have only just re¬ 
cently started investigating coherences in the context of 
resource theories, and a complete understanding is still 
elusive [36]. As we show using the stochastic absorption 


refrigerator, weakly coupled three-qubit fridges operat¬ 
ing in the steady state are in some sense equivalent to 
classical ones, as one can achieve identical steady-state 
temperatures from a classical stochastic model. On the 
other hand, the coherent transport of energy inducing os¬ 
cillations in the population of the cold qubits constitutes 
a unique quantum feature. This points towards the po¬ 
tential of harnessing genuine quantum resources to take 
thermal machines beyond the classically possible. 
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Appendix A: Thermalisation models 

In this Appendix we describe the two models of ther¬ 
malisation used throughout this work, and derive the ap¬ 
propriate master equations. 


1. Model I 

In this model, each heat bath is represented as a col¬ 
lection of harmonic oscillators, so that the total Hamil- 
tonian for the three baths is with 

k 

where the bosonic mode operators satisfy canonical com¬ 
mutation relations ~ ^ij^kk' and = 

0. The qubit-bath interaction is given by 

3 

HAB = Y,Ai®Xi (A2) 

i=l 

where Ai = af and the collective bath coordinates are 
defined by 


k 

with constants A^^k that control the strength of the cou¬ 
pling of qubit i to its associated bath. 


a. Strong-coupling limit 

We now sketch the derivation of the strong-coupling 
master equation, valid when g > Ei. This master equa¬ 
tion describes dissipation as resulting from incoherent 
transitions between the eigenstates of the full coupled 
Hamiltonian Ha = i^ioc + V. When ^ = 0, these eigen¬ 
states are simply the computational basis states. When 
g ^ 0, the interaction splits the degenerate states span¬ 
ning the transport subspace into two new eigenstates, 
denoted by |±) = (|101) ± |010))/\/2, with correspond¬ 
ing energy eigenvalues E 2 Eg. The remaining eigenstates 
and eigenvalues are left unchanged. 

Working in an interaction picture with respect to Ha + 
HB , the time evolution of the system coupling operators 
is given by 


Ai{t) = 

We decompose this into Fourier components as 

[i?^, Ai(w)] = -ujAi{Lo), 

UJ 

(A4) 


where the Bohr frequencies {cj} denote the set of all pos¬ 
sible (positive and negative) energy differences between 
the eigenvalues of Ha^ 


We assume that the initial state of the system 
factorises as /)(0 ) = Pa(0 ) where pB^ = 

/Tr ), and that the system-bath cou¬ 

pling is sufficiently weak that perturbation theory can 
be applied. The master equation is derived by project¬ 
ing the interaction-picture von Neumann equation for the 
density operator onto the subspace spanned by states of 
the form pA ^ pB^ where pA = Tr^(p), and truncating 
the resulting equation at second order in the system-bath 
coupling (Born approximation). The Markov approxima¬ 
tion then consists of assuming that the memory time of 
the bath is much shorter than all typical time scales of the 
reduced qubit evolution, see Ref. [49] for details. Tracing 
over the baths results in the following equation of mo¬ 
tion for the qubit density operator pA in the interaction 
picture: 


dpA 

dt 




UJ,UJ' 2 = 1 


Ai{oj)PA{t)Al{uj') 


- Al{u')Ai{u!)pA{t) 


+ h.c., 


(A5) 


where we defined the self-energy 

poo 

Ti{uj)= dte'‘^^xl{t)Xi{0)), (A6) 

Jo 


with Xi{t) = , In general, the self-energy 

can be written r^(co’) = -h iSi{uj) where the real 

part corresponds to an incoherent transition rate, 

and the imaginary part Si{uj) corresponds to an energy 
shift which we assume to be negligibly small. 


We now perform the rotating-wave approximation by 
averaging over the oscillating terms in Eq. (A5), so that 
terms with uj ^ uj' drop out. This approximation is 
valid when the typical energy differences are much larger 
than the incoherent transition rates, i.e. min{E^,^} ^ 
max{ 7 ^(co’)}. Transforming back to the Schrodinger pic¬ 
ture results in a Lindblad equation of the form 

= —'A^a.Pa] + EE 7 i(w)P[Ai(w)]/)^, (A7) 

2=1 UJ 

where the dissipators are given by 

V[L]p = LpB (A8) 

for a general Lindblad operator L. Explicitly, the Lind- 
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blad operators are 

Ai{Ei) = | 011 ) ( 111 | + | 000 ) ( 100 | 

A(i;i+ 5 ) = ^(|ooi)(+|-|-)(iio|) 
A^{E^-g) = E(^\+) (1101 + 1001 ) (- 1 ) 


wave approximation does not apply to counter-rotating 
terms of frequency {uj — uj') ^ g. Instead we should work 
in an interaction picture generated by H\oc E Hb and 
treat the interaction V between the qubits as a small 
perturbation[30, 31]. The time evolution of the system 
coupling operators is now given by 

Ai(t) = 


AziEz) = | 100 ) ( 110 | + | 001 ) ( 011 | 
A2(E2+g) = ^(i000} (+| + |-)( 111 |) 

A2(E2-g) = ^H+} ( 1111 - 1000 ) (- 1 ) 


^ 3 (-E 3 ) = | 110 )( 111 | + | 000 ) ( 001 | 

A3(E3 + g) = ^(ll00} (+|-|-)( 011 |) 

A3{E3 — g) = 1+) (011| + |100) (—I), 

while the remaining non-zero Lindblad operators, cor¬ 
responding to the reverse processes, are found from 
Ai{-uj) = Ai{ojy. 

In order to actually evaluate the rates we intro¬ 

duce the spectral function of each bath: 

Jj(w) = 27ry]|Ai,k|^<5(w - ^'i.k)- (A9) 

k 

In the limit of an infinite bath with a smooth density of 
states, the sum over the quantum numbers k can be ap¬ 
proximated by an integral, and Ji{uj) becomes a contin¬ 
uous function. We assume that the baths have identical 
spectral functions of the Ohmic form 

J{uj) = (AlO) 


where a is a dimensionless coupling parameter and Q is 
a cut-off frequency of the system-bath interaction, which 
must be much larger than all other energy scales in order 
for the Markov approximation to hold. The incoherent 
rates are then given by 


7i(w) 


J(w)[l + n(w,ft)] (w > 0) .. . 

J{\ioM\u\,l3i) (coKO) ’ 


where n{u,l3) = (e^“ — 1) ^ denotes the Bose-Einstein 
distribution. 


with the corresponding Fourier decomposition 

[Hloc, = -OjAi{oj), 

UJ 


where the frequencies {cj} represent the eigenvalue dif¬ 
ferences of H\oc only. 

As before, we assume that the initial state of the sys¬ 
tem factorises as p(0) = Pa(0) (8) Pb- Now we project the 
interaction-picture von Neumann equation onto states of 
the form pA^ pB and truncate the resulting equation at 
second order in the qubit-bath interaction and the qubit- 
qubit interaction. We then perform the Markov approx¬ 
imation and trace over the bath variables; see Ref. [31] 
for full details of the derivation. The resulting equation 
of motion is 


dpA 

dt 


= -i[V,PA] 

\uj,uj' i=l 


Al{u')Ai{oj)pA{t) +h.c. , 


(A12) 


where the self-energy r^(co’) is defined by Eq. (A6). We 
then write ri(ci;) ^ ^ 7 i(cj), neglecting the imaginary part 
corresponding to small energy shifts of the qubit energy 
splittings. The rotating-wave approximation now con¬ 
sists of crossing off counter-rotating terms with uj ^ uj'^ 
all of which have frequencies of order Ei. In order to be 
consistent with our assumptions, we must have Ei^ g. 
The resulting master equation in the Schrodinger picture 
is 


dpA 

dt 


= - \[Ha,Pa] 


3 


+ E ]pA+li{-Ei)V[af]pA] , 

i=l 

(A13) 


where the rates are given by Eq. (All). 


b. Weak-coupling master equation 


2. Model II 


We now consider the weak-coupling limit where g is 
comparable to the dissipation rates. In this case, the 
previous derivation is no longer valid since the rotating- 


In this model, the baths are represented by three ad¬ 
ditional fictitious qubits described by Pauli operators 
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^ with Hamiltonian 


HF = \'pEiTt. (A14) 

We have chosen these qubits to have identical energy 
splittings to their associated physical qubits, in order to 
avoid renormalising the physical qubit energy splittings. 
We introduce a Lindblad dissipator for each fictitious 
qubit corresponding to damping by a perfectly Marko¬ 
vian (delta-correlated in time) reservoir at temperature 
Ti. For simplicity, we assume identical spontaneous emis¬ 
sion rates 7 for the three fictitious qubits. The coupling 
to the refrigerator is described by the Hamiltonian 

3 

Haf = Ti' + <T~ T+ ), (A15) 

i=l 

where af = ^{crf and = |(rf while 77 is 

a small dimensionless parameter. The density operator 
Paf of the six-qubit system is therefore described by the 
master equation 

= - i[HA + Hf + Haf, Paf] 

3 

+ 7 ^ [Elr^jpAF + e~^*^*'D[TppAF] ■ 

(A16) 

The effective spectral density seen by each physical qubit 
is of the Lorentzian form 




rf + (Ei-c^) 2 ’ 


(A17) 


where Fj = 7(1 + e“^’^‘)/2. For a fixed rj and Tj, the 
bandwidth of the effective bath’s frequency response is 
therefore controlled by modifying the parameter 7 . 


Appendix B: Phase sensitivity analysis 

In this Appendix we consider the effect on cooling 
of coherence in the initial quantum state with a small 
amount of phase noise. In order to gain some analytical 
insight, we suppose that the dynamics can be treated as 
approximately unitary over the time scales of interest. 
This is a decent approximation in the strong-coupling 
limit, where g is much larger than the relaxation rate. 

We define the populations a{t) = (010| |010) 

and b{t) = (101| |101), and the coherence C{t) = 

(1011 PA{t) 1010)5 where pA{t) denotes the reduced quan¬ 
tum state of the refrigerator qubits. We also introduce 
the convenient parametrisations S = [a(0) + 6(0)]/2, 
D = [6(0) — a(0)]/2 and C(0) = ire^"^, where r > 0 is 
the magnitude and (j) is the phase. The refrigeration con¬ 
dition Ty < Ti, where Ty is the virtual temperature (6) 


and Ti is the initial temperature of the cold qubit, implies 
that D > 0. 

At short times, the effect of dissipation can be approx¬ 
imately neglected, and the evolution of the populations 
a(t), b{t) is analogous to that of a resonantly driven two- 
level system. In this approximation, the populations at 
time t are given by 

a{t) = S — D cos(2^t) -h r cos(0) sin(2^t), (Bl) 

and 6 ( 6 ) = 2S — a{t). The population difference of the 
cold qubit is given by 

KW) = K(0)) + 2[a(0)-a(i)], (B2) 

which is related to the effective temperature by Eq. ( 6 ). 
From this we find that the first temperature minimum 
occurs at time 


"“yiv (B3) 

with 

= S+ y/D^ + r^ cos^<P). (B4) 

The amplitude of the population oscillations is therefore 
maximised by choosing 0 = 0 , i.e. a purely imaginary 
initial coherence, as expected from Eq. (7). With the 
choice of 0 = 0 , the optimal time to extract the cold 
qubit is given by topt = 6 min( 0 ). 

We now suppose that each qubit entering the refriger¬ 
ator is prepared with a random phase 0 of the coherence 
in the transport subspace, due to experimental noise, for 
example. These phase shifts are described by some prob¬ 
ability distribution function p(0). In order to find the 
dynamics of the resulting ensemble of cold qubits, we av¬ 
erage Eq. (Bl) over the distribution p(0). This procedure 
leads to the replacement of cos( 0 ) in all expressions by 

cos(0) = f d(/)p(0) cos((/)). (B5) 




As a simple but broadly relevant example, we assume 
that phase shift has zero mean 0 = 0 , and that its vari¬ 
ance 0^ is smaller than unity. We also suppose that its 
cumulants are rapidly decreasing, so that p( 0 ) may be 
approximated by a Gaussian distribution. Under these 
assumptions, we obtain 


cos( 0 ) = e 


= P-^V 2 


(B 6 ) 


We find from Eq. (B3) that the minimum temperature 
of the ensemble occurs at time 


/ 1 / -1 

^opt = ^ 1 ^ 


,-<^ 2/2 

“d 


(B7) 


SO that 6 qp^ > topt in this approximation. In the follow- 













12 


ing, primed variables indicate quantities in the presence 
of noise, while the unprimed variables denote the corre¬ 
sponding quantities in the absence of noise. 

In order to find the shift in the ensemble temperature 
due to the phase noise, we consider two distinct scenar¬ 
ios. In the first scenario, we assume that some informa¬ 
tion about the phase distribution is known, so that it is 
possible to predict in advance the optimal time to 
extract the qubit. This information could be obtained 
by expending some qubits in order to characterise the 
noisy preparation by quantum state tomography. In this 
case, the phase noise leads to a change in the population 
difference of the cold qubit given by 


AK) = K(4J)'-K(topt)) 

= 2Vd2 + r^ -2 \Id^ + r^e-^ 


V-D 2 + r 2 




(B8) 


In the second scenario, we assume that no informa¬ 
tion about the phase noise distribution is known, so that 
the experimenter extracts the cold qubit at the incor¬ 
rect time topt- For example, this scenario applies if only 
a single qubit is available for cooling, so that the noisy 
preparation cannot be characterised in advance. In this 
case, phase noise changes the population difference by 
the amount 


= K(^opt))' - K(^opt)) 


2r^ 

\/D‘^ _1_ ^2 




(B9) 


Note that it is also possible to consider the case where 
only partial information about the noise distribution is 
known, in which case we expect A (erf) to lie somewhere 
between Eqs. (B8) and (B9). 

In any case, assuming that the effect of phase noise 
is small, we find that the expected shift in the ensemble 
temperature to lowest order in A(crf) is 


ATi = ^ cosh^ A) AK), (BIO) 


where Ti denotes the optimal temperature in the absence 
of phase noise. This equation indicates that the effect of 
phase fluctuations is minimised when Ti ^ , while the 

most drastic effects occur for very large or very small 
target temperatures. In particular, the right-hand side 
of Eq. (BIO) diverges as Ti ^ 0, which indicates that 
even vanishingly small fluctuations in the initial coher¬ 
ence phase may prevent one from using this coherence to 
cool the cold qubit to its ground state. 

In Eig. 4(a) we numerically calculate and plot an ex¬ 
ample of the change in the final temperature of the cold 
qubit when it is disconnected from the machine, either 
with or without sufficient knowledge of the phase noise 
distribution to predict the optimal extraction time. We 


(a) (b) 




FIG. 4. Comparison of exact numerical and approximate ana¬ 
lytical results for the effect of Gaussian zero-mean phase noise 
characterised by its variance (jp. The same parameters are 
used as in Fig. 2(a), with 5% of the maximum initial coher¬ 
ence. (a) Shift in the cold qubit temperature when extracted: 
numerical (blue solid line) and analytical (red dashed line) 
results when the phase noise distribution is not known, nu¬ 
merical (black dot-dashed) and analytical (green dotted line) 
results when the phase noise distribution is known, (b) Shift 
in the time of the temperature minimum Atopt = topt ~ Apt: 
numerical (blue solid line) and analytical (red dashed line) 
results. 


see that in fact it makes little difference to the tempera¬ 
ture whether or not the optimal extraction time is known 
precisely, since the shift in this time due to phase noise 
is rather small (Fig. 4(b)). Our approximate analytical 
results are also shown in Fig. 4 for comparison. 


Appendix C: Stochastic refrigerator model 

In this Appendix we derive an effective stochastic 
model of the quantum absorption refrigerator in the deep 
weak-coupling limit. Our analysis follows Ref. [42], where 
it was shown that coherent propagation of a quantum 
particle on a lattice under strong local dephasing can be 
approximated by a stochastic hopping process. In our 
case, the dephasing is provided directly by the action 
of the thermal baths, which destroy the coherences in 
the computational basis at a characteristic rate F. In 
the limit where ^ F, we can derive a closed equation 
of motion for the populations, describing their dynamics 
over times coarse-grained on the scale F“^. 

Working in the weak-coupling limit of Model I, we be¬ 
gin with the master equation (AI3). This may be written 
as 


t = 

where we introduced the Liouvillian superoperator £, 
and for brevity we have written p = pA for the reduced 
state of the refrigerator qubits. We decompose the Li- 
ouvillian as £ = £o + V, where Vp = i[p, E], while the 
superoperator £o contains both the local Hamiltonian 
and dissipative contributions to the Liouvillian. 

The quantum state can be expanded in the eigenbasis 
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of £ as 

p{t) = '^pxe^\ (C2) 

A 

where px is an eigenvector of £ with complex eigenvalue 
A. Since Re(A) < 0, we see that only the eigenvalues with 
the largest real part are relevant as t ^ oc. In this basis 
the master equation reduces to the eigenvalue equation 

(£o + V)pa = Apa. (C3) 

The eigenvalues of jCq are of order T, where T is a char¬ 
acteristic dephasing rate of the thermal dissipation. On 
the other hand, the eigenvalues of V are of order g, which 
is much smaller than T by assumption. We can therefore 
treat the effect of V as a small perturbation. 

We introduce a projector V defined by 

7 

n=0 

where {|n)} are the computational basis states. This 
projects onto the space of populations (diagonal matrix 
elements) in the computational basis. We also define its 
orthogonal complement by Q = 1 — P, which projects 
onto the space of coherences (off-diagonal matrix ele¬ 
ments). We refer to the spaces of populations and co¬ 
herences as the P-space and the Q-space, respectively. 
It is readily verified that the following properties hold: 

VVV = 0, (C5) 

[P,£o] = [Q,£o] =0. (C6) 

Eq. (C5) states that the interaction V only couples pop¬ 
ulations to coherences, while Eq. (C6) reflects the fact 
that the local Liouvillian £o does not couple populations 
and coherences, i.e. it is block-diagonal. 

By introducing the identity 1 = Q-\-V into both sides 
of the eigenvalue equation (C3), we obtain 

XVpx=CoVpx^VVQpx (C7) 

XQpx = ^oQpx + QVQpx + VVpx^ (C8) 

where Eqs. (C5) and (C6) have been used. We now solve 
Eq. (C8) for Qpx and substitute the result back into 


Eq. (C7), finding 

XVpx = CoVpx + VV{\-QV- Co)-^ VVpx- (C9) 

Since we are looking for eigenvalues A which are small 

in magnitude, we may neglect the term A — QV, which 
is of order in comparison to £o, which is of order T. 
Hence we obtain the approximate eigenvalue equation 

XV Px = C,sVpx, (CIO) 

where 

Ceff = Co-VVCpVV. (Cll) 

As t ^ oo, the master equation (Cl) can therefore be 
approximated in the 7^-space by 

^ = £effP. (C12) 

Note that £eff is well defined by Eq. (ClI) despite the 
fact that the unperturbed Liouvillian possesses a zero 
eigenvector (the steady-state solution): jCqPoo = 0. This 
is because the unique steady-state solution lies in the V- 
space, i.e. poo = 'Ppoo- Eurthermore, one can write VV = 
QW by virtue of Eq. (C5). Therefore, the operator 
acts only in the Q-space, where all eigenvalues of £o are 
non-zero. 

In order to calculate £eff explicitly, we note that V 
annihilates all populations apart from |0I0) (0101 and 
|I0I)(I0I|. The action of £efT can therefore be found 
from its action on these two populations. We also make 
use of the fact that 

£o|OIO)(IOI| = -r|0I0)(I0I| 

£o |I0I) (0I0| = -T |I0I) (0I0|, (CI3) 

with 

V = \Y,[lim+ii{-Ei)], (C14) 

i=l 

where ji{Ei) is defined by Eq. (All). Using these prop¬ 
erties, a straightforward calculation yields the result 

CeS = Co + ^V[B] + ^V[B^, (C15) 

where B = |0I0) (I0I|, and V is defined by Eq. (A8). 
This effective Liouvillian describes transport as a 
stochastic process, which transfers population symmetri¬ 
cally between the states |0I0) and |I0I) at a rate 2g^ jV. 


